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Processes mediated by T^-bosons often involve pro- 
duction or annihilation of two fermions with very dif- 
ferent masses. Examples of present interest include the 
_B-meson decay constant fs and the single top-quark pro- 
duction. The basic ingredient in the analysis of such 
processes are correlators of heavy-light currents. For ex- 
ample, Jb can be determined by relating such correlators 
to the measured spectrum of B mesons, with help of the 
QCD sum rules 0. 

In another important application one can use such cor- 
relators, computed perturbatively in the continuum, as 
an input for lattice calculations, in particular for the 
matching of the lattice and continuum currents. The 
perturbative production cross section computed close to 
the threshold (related to the imaginary part of the cur- 
rent correlators) is a convenient physical observable to 
perform the matching B. 




FIG. 1; Examples of diagrams contributing to the current 
correlators at O (aT) . 



Even if the light fermion mass is neglected, a correlator 
still involves two mass scales (see Fig. |l|), the invariant 
mass of the pair i/g^ and the mass of the heavy fermion 
m, which hampers the evaluation of the higher-order 
quantum effects. It is helpful to develop a computational 
scheme which allows to expand the Feynman diagrams 



around the threshold, q 



In the past, various ex- 



pansion schemes were constructed for other kinematical 
situations [§ i |, |, 0, |, §• They facihtated many 
studies of higher-order radiative corrections to a variety 
of processes of experimental importance, such as the Z- 
boson decays, heavy quark production and decays, elec- 
tromagnetic properties of particles, etc. In the present 
paper we present an expansion scheme for the case de- 
picted in Fig. |l|, for (7^ ~ m"^, and apply it to compute 
O (ofg) corrections to the imaginary part of the heavy- 
light correlators of (axial)vector as well as (pseudo)scalar 
currents. The real part can be obtained, if need arises, 
along the same lines. 

Recently, numerical values of those corrections were es- 
timated 0, |ll| using Pade approximants to describe the 
correlator as a function of q^. The approximants were ob- 
tained from several terms in the expansion of the correla- 
tor for large (g^ ^ m^) and small {q^ m^) values of the 
total energy. This approach is accurate for q^ sufficiently 
far from the production threshold q^ = rn? . However, for 
practical applications related to heavy quark physics, one 
needs the correlators quite close to the threshold where, 
for kinematic reasons, their imaginary part vanishes as 
[q^ — m?Y. Because of this suppression, O (a^) terms 
obtained from the Pade approximants are rather uncer- 
tain near the threshold, with an error estimated to be 
about 30% III. 

The approach presented in this Letter enables an ana- 
lytic calculation of those terms. We introduce the param- 
eter 5=1 — rr? jq^ which is small in the kinematic region 
close to threshold. We will show that an expansion in 8 
can be constructed applying the reasoning familiar from 
HQET directly to the Feynman diagrams. The only new 
ingredient necessary for a calculation of the O(a^) cor- 
rections to the correlators are certain three-loop HQET 



2 



diagrams whose calculation we will describe. 
We consider the following correlator, 

(-g'ffM- + q^q.) Tl^'iq^) + q^.q.Ill{q^) 
= 1 [ Axe^'>^0\Tjl{x)fMm, 



(1) 



with = ipij^ip2 and -01,2 denoting the Dirac spinors 
for the massless and massive quarks, respectively. This 
correlator is an analytic function of with a cut starting 
at = m^, corresponding to (5 = 0. Since S is propor- 
tional to the phase space available for the production of 
heavy and light quarks, for small S the heavy quark is 
non-relativistic and always close to its mass shell. On 
the contrary, the massless quark is always ultrarelativis- 
tic, but its energy is small if S is small. The HQET 
is designed to study exactly this kind of situation and 
the calculation of the relevant Feynman diagrams can be 
simplified if one follows its pattern. 

We have to consider two different scales of the momen- 
tum k flowing along a given line in a Feynman diagram: 
hard A: ~ to or soft k ^ mS. Consider the heavy quark 
propagator ^ l/{k'^ + 2qk + q^S). If k is hard, the prop- 
agator can be expanded in a Taylor series in 5 yielding 
the on-shell heavy quark propagator ~ + 2qk). On 

the other hand, if k is soft, the propagator can be ex- 
panded in k^ resulting in the static heavy quark propa- 
gator - l/(2qk + q'^S), familiar from HQET. 

For each diagram, one has to consider all possible mo- 
mentum routings and find all contributing subgraphs. 
Among them, there are two which can be easily de- 
scribed. First there is the situation when all lines are 
soft so that all heavy quark propagators become static 
and the diagram becomes what is usually referred to as 
a HQET matrix element. In our case, these subgraphs 
require three-loop calculations in HQET and we will ex- 
plain below how we solve this problem. 

The second type of subgraphs arises in the situation 
when all momenta are hard. In this case all heavy quark 
propagators are Taylor expanded in i5; the resulting Feyn- 
man diagrams are of the on-shell three-loop propagator 
type, studied in [p^ [l^ , and their evaluation is possi- 
ble. However, since these contributions are polynomials 
in (5, they do not contribute to the imaginary part of the 
correlator and we do not consider them here. 

In between the two extreme cases discussed above, 
there are situations where, in a given diagram, some of 
the lines are soft and some are hard. Using the HQET 
language, these correspond to the HQET matrix elements 
with insertions of higher dimensional operators of the 
HQET Lagrangian or to the HQET matrix elements com- 
puted with leading order operators (HQET currents) cor- 
rected for the higher order Wilson coefficients. Since the 
Wilson coefficients are computed at the hard scale, such 
contributions factorize into products of simple subgraphs 
and can be easily computed. 

Therefore, the main challenge are the three-loop 
HQET diagrams. There are two ways to compute them 



and we have taken both to have a cross check. Recently, 
the three-loop HQET diagrams have been analyzed in 
I p^ , and a computer algebra program has been pub- 
lished, capable of computing all three-loop HQET prop- 
agator type diagrams. We have used that software to 
calculate the required HQET matrix elements. For the 
purpose of the cross check, we have written (in FORM 
|16| ) , in a completely independent way, a similar program 
solving the three-loop recurrence relations for HQET 
propagator-type integrals, restricting ourselves to topolo- 
gies needed for the current calculation. 

In both approaches, every Feynman diagram is ex- 
pressed in terms of a few master integrals. The majority 
of them is known ^ However, one of the master 
integrals has not been evaluated to sufficient accuracy in 
the existing literature and we have to compute it. It turns 
out that one can use a trick to this end. The Euclidean 
integral we need is (p^ = —1, rf = 4 — 2e): 



/ = 



d^fcid^fed^fca 



fc2fc|(fci-fc2)2(fc3-fc2)2 



1 



(2pfci-f I)(2pfc2 + l)(2pfc3 + l)' 
To compute it, consider a similar integral, 



(2) 



d''fcid''/c2d''fc3 



kfk^kliki - k2)Hk3 - k2r 
1 

{2pki + l)i2pk2 + l)i2pk3 + l)' 



(3) 



related to / by integration-by-parts identities. We 
perform a transformation \ki\ — s- l/|fci| for i — 1,2,3 in 
Eq. (0). The integral Ii transforms to 



h = 



d'^fcid'^fc2d''fc3 



(2pfci + kf){2pk2 + fc|)(2pfc3 -I- fc|) ■ 



(4) 



In the next step we notice that Ii is finite in four dimen- 
sions, so that the limit — + 4 can be taken; after that Ii 
becomes equal to one of the on-shell three-loop master 
integrals computed in [Ol. We therefore find: 



Ii^C{e) [271^3 - 5C5] 



(5) 



with C(e) = [tt^ '^r(l + e)]'^. We now use recurrence 
relations to obtain / from Ii and derive: 
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.(6) 



With this integral at hand, all the three-loop HQET mas- 
ter integrals are available to sufficiently high power in 
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FIG. 2: A*'c(5^AX,NA,L,H plotted as a function of = 5/{2-S). 



their expansion in e and we can proceed with the com- 
putation of the correlators. 

We introduce the dimensionless quantity 



T{W* ^ bu) 
\Vbu\^r{W* eP) 



(7) 

and expand it in a series in the strong coupling constant 
(we use the MS scheme and denote a = as{M)/Ti) 



(8) 



We find (an exact formula for R^^^'"" = NcR\ can be 
found in eq. (9)): 
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Rl = (CfA^^ + CaA^V^ + T^TVlA^ + TflA^) (9) 
For the individual color structures we obtain: 
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We denoted Ls = InS and the zeta function ^3 ~ 1.202. 
We do not display higher order terms in the expansion 
in S, but they can be easily obtained as well. 

For completeness, we also give the results for 0{al) 
correction to the correlator of two scalar currents. Let us 
define 

q^U%q') = z J dxe^'^^0\Tf{x)fiOm. (H) 

where js ~ ^mV'iV'2 and, again, '^\;2 denote the Dirac 
spinors for the massless and massive quarks and is the 
on-shell mass renormalization constant for the massive 
quark. We define 



R\(^) 87rlm[ff(q2+i(5)] , 

~ N^{Rl + CFaR\ + CFC?R\ 
We then find (for an exact formula for R\ see [|ll| 



(12) 



Rl = 8\ 
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For the individual color structures we obtain: 
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The results for the pseudo-scalar and axial-vector cur- 
rents are the same as for the scalar and vector currents 
because the presence of massless fermion line permits one 
to cancel the Dirac 75 matrices. 

It is interestin g to compare our results with the num- 
bers obtained in |1^ . We have plotted our results for the 
independent color structures as functions of the velocity 
V = S / {2 — S) using our results for R2, including terms 
O ((5^) which we do not display in (p^). Comparing these 
results with Fig. 4 of ijllj, we find very good agreement 
for V < 0.6 [|l9|; for larger values of v our truncated se- 
ries is not accurate (as can be expected from the very 
nature of expansion) and more terms in the expansion 
are needed. We have also verified the relations between 
the scalar and vector correlators given in eqs. (32-33) in 
0. _ 

In ||ll[, the values of the 0{6^\n 6) terms were esti- 
mated by fitting numerical solutions. Our formulas pro- 
vide analytic results for these coefficients. In the notation 
of eq. (45) of 0, wc find 
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These results agree well with the estimates of [0 for 
the the abehan cff = 21(6) and the fight quark cfl = 
—2.3(7) contributions. The non-abelian part found in 
p4| , Cfa — 1-2(4), differs by about 9 sigma from our re- 
sult, 4.894. Since the production threshold is the most 
difficult place for the Pade approximants it is very likely 
pof that the accuracy of this particular result was over- 
estimated in [O. 



The expansion around the heavy-light threshold pre- 
sented here extends the class of Feynman diagrams which 
can be evaluated analytically. Our approach can be sum- 
marized as applying the HQET directly to Feynman di- 
agrams. As an example application, we computed the 
imaginary part of flavor off-diagonal current correlators 
to 0{aj.) useful as an input for both fs determination 
from QCD sum rules and also for matching of the lat- 
tice and continuum currents. Similar techniques can be 
used to study second order QCD corrections to differen- 
tial distributions in heavy to light semileptonic decays. 
Work on this is in progress. 
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